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Abstract 

It is proved that quasi-exactly soluble potentials (QESPs) corresponding to an oscil- 
lator with harmonic, quartic and sextic terms, for which the n + 1 lowest levels of a given 
parity can be determined exactly, may be approximated by WKB equivalent potentials 
corresponding to deformed anharmonic oscillators of SU 9 (1,1) symmetry, which have been 
used for the description of vibrational spectra of diatomic molecules. This connection 
allows for the immediate approximate determination of all levels of the QESPs. 
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Quantum algebras [1,2] (also called quantum groups), which are deformed versions of 
the usual Lie algebras, to which they reduce when the deformation parameter is set equal 
to unity, have recently been attracting considerable attention. The interest in possible 
physical applications was triggered by the introduction of the (/-deformed harmonic oscil- 
lator by Biedenharn [3] and Macfarlane [4] in 1989, although similar oscillators had already 
been in existence [5-7] . The common mathematical structure of these various types of os- 
cillators was shown by several authors (see [8] for a list of references), using the concept 
of the generalized deformed oscillator [9] . 

One way to clarify the physical content of deformed oscillators is to try to construct 
potentials which give, exactly or approximately, the same spectrum as these oscillators. 
Along these lines, WKB equivalent potentials corresponding to the (/-deformed oscillator 
[10], as well as to deformed anharmonic oscillators [11] found useful in the description of 
vibrational spectra of diatomic molecules [12,13] have been constructed, both numerically 
[10] and in analytic form [11]. These potentials give approximately the same spectrum as 
the deformed oscillators under discussion. 

On the other hand, several quasi-exactly soluble potentials (QESPs) have been intro- 
duced [14-18], for which the exact calculation of the first n + 1 energy levels is possible, 
with no information provided for the rest of the levels. In particular, in the case of the 
harmonic oscillator with both quartic and sextic anharmonicities [14] , the exact calculation 
of the first n + 1 energy levels of a given parity is possible. 

It is therefore interesting to check if there is any relation between the WKB equivalent 
potentials (WKB-EPs) approximating the behavior of the deformed oscillators on the one 
hand and the quasi-exactly soluble potentials on the other. If such a link exists, as we shall 
demonstrate in this letter, one can fix the parameters of a WKB potential equivalent to 
an appropriate oscillator so that the potential agrees, up to the sextic term, with a given 
quasi-exactly soluble potential. If the approximations involved (WKB method, truncation 
of the WKB potential to sixth order only) are good, one can subsequently use the levels 
of the relevant oscillator in order to approximate the unknown higher levels of the quasi- 
exactly soluble potential. In the case of the harmonic oscillator with both quartic and 
sextic anharmonicities [14], the levels of the opposite parity can also be approximated in 
this way. 

The q-numbers are defined as 



In the case that q = e T , where r is real, they can be written as 

M = (2) 
sinh(r) 

while in the case that q is a phase (q = e lT , with r real) they take the form 

W = (3) 
sm(r) 

It is clear that in both cases the q-numbers reduce to the usual numbers as q — > 1 (r — > 0) . 

The q-deformed harmonic oscillator [3,4] is determined by the creation and annihila- 
tion operators a + and a, which satisfy the relations: 

aa + -q Tl a + a = q ±N , [N,a + ]=a + , [N, a] = —a, (4) 

where N is the number operator. The Hamiltonian of the q-deformed oscillator is 

H = —(aa + + a + a), (5) 



and its eigenvalues are [3,4] 



E(n) = — ([n] + [n+l]). (6) 



From Eq. (6) it is clear that the spectrum of the q-deformed harmonic oscillator is 
not equidistant. For real q the spectrum increases more rapidly than in the classical case, 
while for q complex it increases less rapidly than in the classical case, i.e. it gets squeezed. 

Using two q-deformed oscillators described by the operators ai, a\ , Ni, and 02, aj, N2, 
the generators of the quantum algebra SU g (l,l) are written as [19,20] 

K+ = a+a+, K =ai a 2 , K = \ (N x + N 2 + 1) , (7) 

and satisfy the commutation relations 

[K , K±] = ±K ± , [K+, K_] = - [2K ] , (8) 

where the square bracket in the rhs of the last equation is a q-number as defined in eq. 
(1). In ref. [12] this symmetry has been used for the description of vibrational molecular 
spectra. The q-deformed anharmonic oscillator used there is a q-generalization of the 
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anharmonic oscillators used in the usual Lie algebraic approach to molecular spectroscopy. 
(See ref. [12] for detailed references.) 

The energy spectrum of the anharmonic oscillator with SU g (l,l) symmetry, in the 
case of complex q, is given by [12] 



/t sin(r(n-iV/2))sin(r(n + l-iV/2)) 

b n = b Q -A . 2/ x ' yV) 

sin (t) 



where 

N = 2n max or N = 2n max + 1, (10) 

with n m ax corresponding to the last level before the dissociation limit. For r^Oa Morse 
[21] or modified Poschl-Teller [22,23] spectrum is obtained. The WKB equivalent potential 
in this case is [11] 

-r r / V -r r A / T Si U ( NT ) \ 2 J, 

V( X ) = V m in + - Ay '- U 2 
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sin 2 r 



_ 2 r 2 cos(A^r) 2 1 r 4 (23 cos 2 (iVr) - 6) 



.4 



3 sin 2 r 45 sin 4 r 

u 6 + 



2 r 6 (67 cos 2 (Nr) - 36) cos( A^r ) ^ , ; 



315 ~;- 6 



sin r 



(11) 



where 



and 



V min = E' -A m ^l:T (TN \ d2) 
2 sin r 



u = — - — x. (13) 

It can be easily seen that for r — > Eq. (11) reduces to the Taylor expansion of the 
modified Poschl- Teller potential 



V(x) = V min + ju 2 



1 2 2 17 4 62 6 

1 u 2 H w 4 w b + 

3 45 315 



(14) 



where u is given by eq. (13), while the modified Poschl-Teller potential [22,23] in closed 
form is 

V(x) = V min + — tanh 2 J . (15) 

This means that the WKB-EP of a system with SU q (l, 1) symmetry is a deformation of 
the modified Poschl-Teller potential. 
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On the other hand, it is known [14,15] that the potentials 



V(x) = 8aV + 8abx 4 + 2[b 2 - (2k + 3)a]x 2 , (16) 

with k = 2n + r, where n = 0, 1, 2, . . . and r = 0, 1 are quasi-exactly soluble. The meaning 
of this term is the following: for these potentials one can construct exactly the first n + 1 
levels with parity (— l) r . (The extra factors of 2 in eq. (16) in comparison to refs [14,15] 
are due to the fact that we use the usual form of the Schrodinger operator 



while in refs [14,15] the form 



is used.) 

We wish to check if the WKB-EP for the deformed oscillator mentioned above corre- 
sponds, up to terms of x 6 , to quasi-exactly soluble potentials of the form of eq. (16). For 
reasons of convenience in eq. (13) we set % = m = 1. Comparing the coefficients of x 2 , x 4 
and x 6 in eqs (11) and (16) we get respectively the following equations 

/ 9 / , . . A 2 t 2 sin 2 ( Nt) , , 

2 (6 2 - (2k + 3 a = — r-t (19) 

2 sin t 



3 sin t 



8a 2 = 2 A , r Q sin 2 (Nr) (23 cos 2 (Nr) - 6) _ ( 

45 sin 8 r 

We can consider a 2 = 1 without loss of generality. For the ground state to be normalizable 

one should have a > [16], which implies in the present case a = 1. Having in mind that 

in eq. (9) one must have A > in order to get an increasing spectrum (this fact can be 

seen also in the realistic applications to molecular spectra given in ref. [12]), from eq. (21) 

one finds 

A = JeV^ - - S in / (22) 

V r 3 / 2 v /sm(NT){/23 cos 2 (Nr) - 6 

Because of the symmetry q <-> g _1 characterizing the q-deformed oscillators, it suffices to 
consider r > 0. The following conditions should then be satisfied 

sin(iVT) > 0, (23) 
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cos 2 (iVr) > 6/23. 



(24) 



Then eq. (20) gives 



b 



15y/E A/sin(A^r) cos(iVr) 



(25) 



2 V?(23cos 2 (iVT)-6) 3 / 4, 



while eq. (19) gives 



2k + 3 = 



3^ sin(A^r) (18 cos 2 (Nt) - 6) 
~2 r(23cos 2 (A^r) -6) 3 / 2 



(26) 



The fact that k = 2n + r has to be positive gives the additional condition 



cos 2 (iVT) < 1/3. 



(27) 



From eq. (25) it is clear that b can be either positive (when cos(iW) < 0) or negative 
(when cos(iW) > 0). 

For reasons of convenience, let us consider the case b > first. In this case the 
conditions sin(iW) > and cos(iW) < imply that we can limit ourselves to the region 
7T > Nt > tt/2. The conditions (24) and (27) then imply that 2.1863 + 2nl > Nt > 
2.1069 + 2nl, with / = 0, 1, 2, . . . . Thus though we can find an infinity of values for N, as 
we shall see it suffices to consider / = 0. 

We first wish to check the accuracy of our approach, which already involves two 
major approximations: the WKB approximation and the omission in the WKB-EP of 
terms higher than x 6 . In order to achieve this, we will compare the results given by this 
method to the exact results given in ref. [14] for the case n = 1, r = 0, in which 2k + 3 = 7. 
From eq. (26) one then sees that a possible solution is iV = 151, r = 0.0144503, which gives 
A = 0.4343473, b = 12.589097, while from eq. (12) (for V min = 0) one has E' = 1636.8943. 
We see therefore that eqs (19)-(21) are indeed satisfied and the potentials of eqs (11) and 
(16) (up to the x 6 term) are identical, being 



E = 3b-2Vb 2 + 2 = 12.4307, E 2 = Sb+2y/tf + 2 = 63.1039, E 2 -E = 50.6731, 

(29) 



V(x) = 303.02(x 2 + 0.3324x 4 + 0.02640x 6 ). 



(28) 



The two solutions which can be obtained exactly are [14] 
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while eq. (9) gives the complete spectrum, including 

E = 12.3805, E 2 = 63.0584, E 2 - E = 50.6779. (30) 

We remark that the agreement between the exact results of ref. [14] and the predictions 
of eq. (9) is excellent, implying that in this case: 

i) the WKB method is accurate, 

ii) the omission of terms higher than x G is a good approximation. 

An additional check for the case n = 3, r = is also made. In this case 2k + 3 = 15. 
One possible solution of eqs (19)-(21) in this case is given hy N = 325, r = 0.00671384, 
A = 0.2960795, b = 18.469158, E' Q = 5168.941. The potential is 

V(x) = 652.20(x 2 + 0.2265x 4 + 0.01227x 6 ). (31) 

Eq. (9) gives the levels 

E = 18.1126, E 2 = 91.3482, E 4 = 165.8771, E G = 241.6456, (32) 

while, following the procedure of ref. [14] , one sees that in this case the exact energy levels 
are the roots of the equation 

E 4 - 28bE 3 + (2546 2 - 240)£ 2 + (-8126 3 + 25926)£ + (5856 4 - 46566 2 + 2880) = 0, (33) 

given by 

E = 18.1429, E 2 = 91.3783, E 4 = 165.913, E 6 = 241.703. (34) 

We remark that excellent agreement between the approximate and the exact results is 
again obtained. 

The agreement remains equally good if more levels are considered. For the case n = 9, 
r = 0, in which 2k + 3 = 39, one possible solution of eqs (19)-(21) is given by iV = 399, 
T = 0.00545864, A = 0.2703882, b = 21.275801, E' Q = 7126.0336. The potential is 

V(x) = 827.43(x 2 + 0.2057x 4 + 0.009669x 6 ). 

The energy levels are compared to the exact ones in Table I. 
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We turn now to the case with b < 0. The conditions sin(A^r) > and cos(iW) > 
imply that we should limit ourselves to the region 7r/2 > Nt > 0. Then the conditions 
(24) and (27) imply that 1.0347 > Nt > 0.9553. (We again ignore terms of 2nl.) 

We again consider the case n = 1, r = 0, 2k + 3 = 7. One possible solution is given by 
N = 61, r = 0.0157377, A = 0.4538508, b = -12.108743, E' = 390.66689. The relevant 
potential is 

V(x) = 279.095(x 2 - 0.3471x 4 + 0.02866a: 6 ). (35) 
The exact solutions given by ref. [14] are 

E = -60.7083, E 2 = -11.9441, (36) 

while eq. (9) gives 

E = 11.7456, E 2 = 57.3764. (37) 

The reasons for this dramatic failure are well understood. The potential of eq. (35) has 
a central well in the middle, having its minimum at x = 0, E = 0, plus two symmetrically 
located wells, on either side of the central one, with their minima at E < 0. The method 
of ref. [14] gives the lowest two energy eigenvalues of the system, which in this case are the 
lowest two levels in the side wells. The WKB method used in ref. [11], however, is known 
to be valid only for small values of x, i.e. in the area of the central well. Therefore in this 
case eq. (37) gives the lowest two levels of the central well, which in this case are not the 
lowest two levels of the system. 

We therefore conclude that a correspondence between the SU 9 (1,1) WKB-EP and the 
quasi-exactly soluble potentials of eq. (16) can be made only for b > 0, in which case both 
methods will give the levels corresponding to a well around x = 0. In this case, as we have 
already seen, the approximation is very good, providing us with the following method: 

Given a quasi-exactly soluble potential for which only the n + 1 lowest levels of parity 
(— l) r can be found exactly, we can choose appropriately the parameters of an SU g (l,l) 
WKB-EP, so that the two potentials are identical up to order x 6 . Then using eq. (9), i.e. 
the known eigenvalues of the WKB-EP, we can approximate the levels of the same parity 
lying above the first n + 1 known levels, as well as the levels of opposite parity, for which 
the method of ref. [14] gives no information for the potential under discussion. (It is worth 
noting that changing the values of n and r changes the potential.) 

Using a similar procedure one can easily see that no connection is possible between 
the WKB-EPs [11] of the q-deformed harmonic oscillator and the QESPs. In the case of 
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the Q-deformed harmonic oscillator [5-7], the WKB-EP of which has been given in [24], 
no connection is possible either. The same is true for the modified Poschl- Teller potential, 
the Taylor expansion of which is given by eq. (14). 

The following comments on the results can now be made: 

i) The quasi-exactly soluble potentials are known to be related to the SL(2) [15] 
and deformed SU(2) [25] symmetries. Among the several oscillators mentioned here, the 
only one for which the connection to QESPs was possible is the oscillator with SU 9 (1,1) 
symmetry. Apparently the H 9 (4) symmetry of the q-deformed oscillator is "not enough" 
to guarantee such a connection. 

ii) The QESPs of eq. (16) are characterized by 3 parameters (a, 6, k). The SU g (l,l) 
potential is also characterized by 3 parameters (A, r, N), while the WKB-EPs of the q- 
oscillator and Q-oscillator are characterized by 2 parameters (u>, r and uj, Q respectively) 
and the modified Poschl- Teller potential is characterized by only one parameter (A). It is 
therefore not surprising that a connection is possible only for the case in which the number 
of parameters of the two potentials to be related is the same. 

In summary, we have proved that the quasi-exactly soluble potentials corresponding 
to a harmonic oscillator with quartic and sextic anharmonicities can be approximated by 
the WKB equivalent potentials corresponding to a deformed anharmonic oscillator with 
SU g (l,l) symmetry. As a result one can use the relevant SU g (l,l) oscillator in order to 
approximate the levels of the QESP which cannot be obtained using that approach. In 
the specific case under discussion the extra levels obtained in this way are the levels of the 
same parity lying above the first n + 1 known levels, as well as all the levels of the opposite 
parity. The fact that the modified Poschl- Teller potential with SU(1,1) symmetry cannot 
be connected to these QESPs, while its deformed version with SU g (l,l) symmetry can, is 
an example where the use of q-deformation allows the solution of an otherwise intractable 
problem. 

In ref. [12] the deformed anharmonic oscillator with SU g (l,l) symmetry has been 
proved appropriate for the description of vibrational spectra of diatomic molecules. An 
extension to vibrational spectra of highly symmetric polyatomic molecules has been given 
in [26]. The construction of QESPs related to realistic vibrational spectra is receiving 
attention. 

One of the authors (DB) has been supported by CEC under contract ERBCHBGCT 
930467. Another author (HM) acknowledges support by KFUPM and the warm hospitality 
during his stay in NCSR "Demokritos" . 
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Table 1 

Exact energy levels of the quasi-exactly soluble potential of eq. (16) with n = 9, r = 0, 
a = 1, b = 21.275801, compared to the approximate energy levels of the WKB-EP potential 
of eq. (11) with A = 0.2703882, E' Q = 7126.0336, r = 0.00545864, N = 399, given by eq. 
(9). 

n E n (approx) E n (exact) 






20.3991 


20.4153 


2 


102.672 


102.682 


4 


186.131 


186.138 


6 


270.735 


270.749 


8 


356.444 


356.485 


10 


443.217 


443.317 


12 


531.013 


531.218 


14 


619.791 


620.165 


16 


709.507 


710.133 


18 


800.118 


801.101 
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